Wavepacket scattering on graphene edges in the presence of a (pseudo) magnetic field 



D. R. da Costa, 1 A. Chaves, 1 G. A. Farias, 1 L. Covaci, 2 and F. M. Peeters 2 ' 1 

1 Universidade Federal do Ceard, Departamento de Fisica Caixa Postal 6030, 60455-760 Fortaleza, Ceard, Brazil 
2 Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium 

The scattering of a Gaussian wavepacket in armchair and zigzag graphene edges is theoretically 
investigated by numerically solving the time dependent Schrodinger equation for the tight-binding 
model Hamiltonian. Our theory allows to investigate scattering in reciprocal space, and depending 
on the type of graphene edge we observe scattering within the same valley, or between different 
valleys. In the presence of an external magnetic field, the well know skipping orbits are observed. 
However, our results demonstrate that in the case of a pseudo-magnetic field, induced by non-uniform 
strain, the scattering by an armchair edge results in a non-propagating edge state. 
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I. INTRODUCTION 

Due to its unique electronic properties, graphene has 
become a topic of intensive study in recent years. Within 
the low energy approximation for the tight-binding 
Hamiltonian of graphene, electrons behave as massless 
Dirac fermions, with a linear energy dispersion.^ This 
leads to a plethora of interesting physical phenomena, 
ranging from Klein tunneling and other quasi-relativistic 
effects 2 3 to the existence of new types of electron de- 
grees of freedom, namely, the pseudo-spin, related to the 
distribution of the wave function over the carbon atoms 
belonging to the different triangular sub-lattices compos- 
ing the graphene hexagonal lattice, and the presence of 
two inequivalent electronic valleys, usually labeled as K 
and K', in the vicinity of the gapless points of the energy 
spectrum of graphene. 

Recent papers studied the scattering of electrons by 
edges 4 and defects 5 in graphene, both theoretically 6 
and experiment alljPl Armchair and zigzag are the 
two types of edges which are most frequently consid- 
ered in the study of graphene ribbons, although other 
types of terminations exist due to edge reconstruction, 
which has been demonstrated both theoretically and 
experiment alf^El. Even so, the edge reconstruction ef- 
fect strongly depends on how the nanoribbon is made: 
normally, it occurs when the technique used to fabri- 
cate the nanoribbon is based on a mechanism that drives 
the system to thermodynamic equilibrium. According 
to the continuum (Dirac) model, armchair edges in fi- 
nite graphene samples lead to a boundary condition that 
mixes the wavefunctions of K and K' valleys, whereas 
a zigzag edge appears in the Dirac theory of graphene 
as a separate bou ndary condition for the wavefunctions 
of each valleyJ^UHl This suggests that electrons reflected 
by a graphene edge would exhibit inter-valley scatter- 
ing only in the armchair case, whereas reflection by a 
zigzag edge would produce scattering inside the same 
Dirac valley. This prediction was confirmed by recent 
experiments where inter- valley scattering by armchair 
edges was even shown to be very robust in the presence 
of defects. The inter- and intra- valley scattering possi- 



bilities are schematically illustrated in Fig. 1(a), which 
shows K and K' Dirac cones in the reciprocal space of 
graphene. 

Besides its singular electronic properties, graphene 
also exhibits interesting mechanical properties, as it can 
support strong elastic stretch. This provides us with 
the new possibility to tune the electron properties in 
graphene through strain engineering In fact, it has 
been demonstrated recently that electrons in a strained 
graphene lattice behave as if they were under an exter- 
nal magnetic field, which points towards opposite direc- 
tions in the K and K' valleys, so that the time reversal 
symmetry of the system as a whole is preserved! 2 ^ Such 
fields were experimentally observed recently, when mea- 
surements of the energy states in a graphene bubble re- 
vealed a Landau level-like structure corresponding to an 
external magnetic field of ~ 300 Tp3 By designing non- 
uniform strain fields in a graphene sheet, one is able to 
produce a uniform pseudo-magnetic field for electrons.^ 
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FIG. 1: (Color online) (a) Dirac cones of graphene, along with 
an illustrative scheme of the inter- (green circles) and intra- 
(gray circles) valley scattering, (b) Sketch of the strained 
graphene sample considered in this work, where the (open 
green) full black circles represent the (un)strained case. The 
upper boundary is set as the y = axis for convenience. 

The aim of this paper is two fold - we use wavepacket 
dynamics calculations: (i) to investigate electron reflec- 
tion by armchair and zigzag edges in a finite graphene 
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sample, assumed to be made by cutting a graphene mono- 
layer, such that no edge reconstruction is expected to oc- 
cur at room temperature, where our results demonstrate 
the possibilities of inter- and intra-valley scattering, de- 
pending on the type of edge, and (ii) to study the in- 
fluence of an external magnetic field and a non-uniform 
strain distribution on the electron trajectories in these 
systems. We compare the features observed for electrons 
under a perpendicular external magnetic field with those 
seen with a pseudo-magnetic field. Fig. 1(b) shows a 
sketch of the graphene flake considered in our calcula- 
tions, where the open (green) circles illustrate the un- 
strained sample and the closed (black) circles illustrate 
the strained one. Such a non-uniform strain field was 
suggested by Guinea et alW* and was shown to exhibit 
an almost uniform pseudo-magnetic field. 

All the calculations were done within the tight-binding 
description of graphene, using the time-evolution method 
developed in Ref. [29J As we are not restricting our- 
selves to a single Dirac cone in our model, the scatter- 
ing between Dirac cones by armchair edges will appear 
naturally. Notice that Fig. 1(b) is just an illustrative 
scheme of our system, where the number of atoms was 
reduced in order to help its visualization. Besides, the 
sample shown in Fig. 1(b) is a ribbon, which improves 
the visualization of the strained case. However, the ac- 
tual flake considered in our calculations has 1801x2000 
atoms, which looks more like a rectangle, rather than 
a ribbon and corresponds to a flake with dimensions of 
about 426x221 nm 2 . Such a large flake is necessary to 
isolate each reflection of the wavepacket on a single edge, 
as we nee d to co nsider a large packet in order to avoid 
dispersion™Sfl 



II. EDGE-DEPENDENT SCATTERING 

Let us first analyze the wavepacket reflection by zigzag 
and armchair edges in a plain graphene sample, i.e. in 
the absence of magnetic fields and strain. The initial 
wavepacket ^(x,y) is taken as a circularly symmetric 
gaussian distribution of width multiplied by a plane 
wave with wave vector k — (k x , k y ) and a pseudo-spinor 
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The pseudo-spinor in our model is simulated by defin- 
ing a multiplication factor in the wave function, which 
assumes different values for sites belonging to the A and 
B sub-lattices. Notice that for low energy electrons in 
graphene, the tight-binding Hamiltonian can be approx- 
imated as H = vpfih • <t, where vp is the Fermi velocity, 
so that the propagation velocity vector in the Heisenberg 
picture is given by dx/dt = —[x,H]i/h = vpS. Hence, 



the pseudo-spin polarization of the wavepacket plays an 
important role in defining the direction of propagation. 
As the upper and right edges of the flake are of armchair 
and zigzag type, respectively, we consider a = (l,z) T , 
i.e. propagation in the ^/-direction, in order to observe 
wavepacket scattering on the upper armchair edge, and 
a = (1, 1) T , i.e. propagation in the x-direction, for scat- 
tering on the right zigzag edge. The initial wave vec- 
tor k is taken in the vicinity of the Dirac point K = 
(0, 47t/3a/3&), where a = 1.42 A is the inter-atomic dis- 
tance. At each time step, we calculate the average values 
( x ) = I-oo x \^\ 2 dxdy and (y) = y\^\ 2 dxdy, in order 
to track the wavepacket trajectory in real space. Besides, 
a fast Fourier transform (FFT) of the wavepacket is taken 
at each time step, in order to track its scattering in re- 
ciprocal space. 

The average positions (x) (black solid) and (y) (red 
dashed) are shown in Fig. 2 as a function of time, for a 
wavepacket propagating in the x (y) direction, towards 
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FIG. 2: (Color online) Average position of a gaussian 
wavepacket of width d = 300 A in an unstrained graphene 
flake, in the absence of external magnetic fields, as a function 
of time, (a) Horizontal propagation of a wavepacket with 
k = (0.03A _1 , 47r/3V3a), and its consequent reflection by a 
zigzag edge. In this case the origin of the system is shifted, 
so that the right zigzag edge is set as the x = axis. The 
wavepacket starts at (xo,2/o) = (— 600A , — 600A ), and ex- 
hibits also a slow drag towards the upper edge, due to Zitter- 
bewegung effects, (b) Vertical propagation of a wavepacket 
with k = (0,47r/3>/3a + 0.02A" 1 ), starting at (0, -600A), 
exhibiting reflection by the upper armchair border. In this 
case, the sample is not shifted, i.e. the upper edge is at the 
y = axis, as sketched in Fig. 1(b). 
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the right zigzag (upper armchair) edge of the sample, and 
being reflected by this edge back to its initial position. 
Figs. 2(a) and (b) correspond to zigzag and armchair 
reflections, respectively. The wavepacket starts at 600 A 
from the sample edge and reaches the edge at t « 70 fs in 
both cases. Due to the finite width of the packet (d = 300 
A ), its center of mass never reaches the border^, so that 
(x) or (y) start to exhibit backscattering when they are 
still ~ 150 A far from the edge. Notice that the motion 
in the y— direction shown in Fig. 2(b) is perfectly verti- 
cal, i.e. (x) = during the whole propagation. However, 
this is not the case for propagation in the x— direction as 
apparent in Fig. 2(a), which is not perfectly horizontal, 
i.e. (y) does not stay the same, as the wavepacket slowly 
drags towards larger y during propagation. This effect 
is a manifestation of the zitterbewegung, as discussed in 
detail in Ref. 29. Although we did not manage to con- 
struct a wavepacket that propagates perfectly horizontal, 
avoiding such a vertical drag, this effect does not interfere 
in our results and conclusions, as our analysis of scatter- 
ing on the zigzag edge depends only on the horizontal 
component of motion. 

Once we know the instant when the wavepacket is re- 
flected by the graphene edge in real space, at that mo- 
ment we analyze what happens in reciprocal space. Fig. 
3(a) shows the lines (red dashed) in reciprocal space along 
which we will take the wave functions. The contour plots 
in Fig. 3(b) illustrate the wave function along the hor- 
izontal line (i) — (ii) depicted in Fig. 3(a) in reciprocal 
space, as time elapses, in the case of ^-direction propaga- 
tion and, consequently, zigzag edge reflection. For such a 
propagation direction, we assumed the initial wave vec- 
tor as k = (0.03A _1 , 47r/3v // 3«). Therefore, the initial 
wavepacket (at t = 0) has a peak around k l ~ l% = 0.03 
A _1 . This peak is conserved until the wavepacket starts 
to be reflected by the right zigzag edge, when interfer- 
ence patterns start to show up. At « 70 fs, a peak 
at k % ~ %% = —0.03 A -1 starts to appear, while the for- 
mer peak at k % ~ %% = 0.03 A -1 smoothly decays. This 
is indeed the instant when the wavepacket is reflected 
by the zigzag edge in real space, as shown in Fig. 2(a). 
As time elapses, the wavepacket ends up only around 
k % ~ %% = —0.03 A _1 . This is direct evidence of intra- valley 
scattering as schematically illustrated in Fig. 1(a). 

Figs. 3(c) and (d) show the wave function in reciprocal 
space taken along the vertical line (Hi) — (iv) depicted 
in Fig. 3(a) as time elapses, for vertical propagation 
and, consequently, armchair edge reflection. For propa- 
gation in the ^/-direction, we consider k = (0, 47r/3v // 3« + 
0.02A -1 ), so that the wavepacket initially exhibits a 
peak around ~ 1.723 A -1 , as shown in Fig. 3(c). This 
peak is preserved up to t ~ 70 fs, when the wavepacket 
is scattered by the upper armchair edge (see Fig. 2 (b)) 
and the amplitude of the peak starts to decrease. Mean- 
while, another peak appears around ky l ~ lv ~ —1.723 
A -1 , which is located in the K' valley, as shown in Fig. 
1(a). The inter- valley scattering situation is illustrated 
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FIG. 3: (Color online) Time evolution of the wavepacket in 
reciprocal space corresponding to the situations shown in Fig. 
2. (a) Illustrative scheme of the lines in reciprocal space along 
which Fourier transform of the wave functions are taken. For 
the propagation in the horizontal direction (see Fig. 2(a)), we 
consider k = (0.03A _1 , 4ir / 3y/3o>) • The time evolution of the 
wave function along the (i) — (ii) line of reciprocal space is 
shown in (b) as contour plots. For vertical propagation (see 
Fig. 2(b)), we consider k = (0, 4vr/3^3a + 0.02A _1 ). The 
time evolution of the wave function along the (Hi) — (iv) line 
of reciprocal space is shown as contour plots in (c) and (d), 
corresponding to different ranges of k % y % ~ %v . 



by the green circles in Fig. 1(a) which is clearly observed 
in reciprocal space. 



III. SKIPPING ORBITS 

Let us now investigate the trajectory of a wavepacket in 
the presence of an external magnetic field, while it under- 
goes reflection at the edges of our rectangular graphene 
flake. We consider the same conditions as in Fig. 2(a), 
i.e. the wavepacket in this case moves to the right, being 
thus pushed to the upper armchair edge by the Lorentz 
force due to the perpendicular magnetic field. The tra- 
jectory drawn by r= ((#), (y)) for such a packet in the 
xy — plane after a t = 2000 fs propagation, under a ^5 
T field, is shown in Fig. 4(a), where skipping orbits are 
clearly observed ^coming from the successive reflections 
at the borders of the system, followed by ciclotronic semi- 
circles, as one would expect from such a scattering prob- 
lem. The arrows indicate the direction of propagation, 
and the edges of the figure are set to be exactly at the 
position of the edges of the graphene flake. In order to 
help their analysis, the trajectories where divided into 
four regions, labeled from (I) to (V). Fig. 4(b) shows (x) 
and (y) separately as a function of time, where one veri- 



4 




(V) 




— d = 100 A 

- - d = 300 A 



-2000 -1000 1000 2000 

x(A) 



s 

o 



o 
a 

- 

0tf 



(IV) ; (V) 

d = 100 A 

- - d = 300 A 




500 1000 1500 2000 

t(fs) 

FIG. 4: (Color online) Trajectories drawn by (x) and (y) for a 
gaussian wavepacket in the presence of an external magnetic 
field w 5 T, propagating close to the edges in a rectangular 
graphene flake, within a t = 2000 fs propagation time. The 
arrows indicate the direction of propagation. Different values 
of the initial wavepacket width d are considered. The edges 
of this panel are placed at the positions of the actual edges 
of the sample, (b) Average values of the wavepacket position 
(x) and (y) as a function of time for the trajectories drawn 
in (a). Different parts of the trajectory in (a) were labeled 
from (I) to (V), and the time intervals where they occur are 
delimited by the vertical lines in (b). 



fies e.g. the attachment of the packet to the upper edge 
((y) close to y = in region (I)), followed by a decrease 
in (y), when it attaches to the right edge ((x) close to 
x = 2100 A, in region (II)), and its further attachment 
to the bottom edge of the sample ((y) close to y = —2100 
A in region (III)). As previously mentioned, due to the 
finiteness of the packet width, the trajectory as described 
by ((x), (y)) does not reach the edges of the system. Be- 
sides, the wavepacket disperses as time elapses, which 
distorts the trajectory as compared to the one obtained 
by classical ballistic motion.^ Even so, the main conclu- 
sion one draws from this result is quite clear: as well as in 
ordinary systems with confined Schrodinger particles,^ 
electrons in graphene under external magnetic fields ex- 
hibit a skipping orbit pattern when propagating close to 
the edges of the sample. We performed calculations for 
different wavepacket widths d = 100 A and 300 A , and 
the results lead to the same qualitative conclusion, dif- 



fering only by the distance the wavepacket may reach the 
edge. There is, however, an important difference between 
these skipping orbits and those in ordinary Schrodinger 
systems, namely, the wavepacket in this case may scatter 
not only between momentum states with opposite signs 
within the same valley, as usual, but they can also scat- 
ter from one valley to another, depending on the type 
of edge, as we demonstrated in Fig. 3. Nevertheless, 
the effect of an external magnetic field on electron states 
belonging to both valleys is the same, therefore, there 
is no detectable manifestation of inter-valley scattering 
in this situation. This is not the case when, instead of 
an external magnetic field, we consider a strain induced 
pseudo-magnetic field, as we will demonstrate in what 
follows. 
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FIG. 5: (Color online) (a) Trajectories drawn by (x) and (y) 
for a t = 2000 fs time evolution of a wavepacket, which prop- 
agates close to the upper (armchair) border of a bent rect- 
angular graphene sample, for two values of wavepacket width 
d. The radius of the circular distortion is R — 10 4 A , corre- 
sponding to an almost uniform w 5 T pseudo-magnetic field. 
The horizontal dashed line represents the upper edge of the 
sample, (b) Average values of the wavepacket position (x) 
and (y) as a function of time for the trajectories drawn in (a). 

For a circularly strained graphene flake, like that 
sketched in Fig. 1(b), electrons in the sample behave as 
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if they were in an almost uniform magnetic field perpen- 
dicular to the plane. In order to produce such a strain, 
our 1801 x 2000 atoms sample is distorted into a semi- 
circle of radius R = 10 4 A , leading to a pseudo-magnetic 
field « 5 T, i.e. close to the value considered for the 
external magnetic field in Fig. 4. The presence of such a 
pseudo-magnetic field when electrons move close to the 
edge are expected to result in skipping orbits, similar to 
those in Fig. 4. Surprisingly, Fig. 5(a) shows this is not 
really the case: after performing a semi-circular trajec- 
tory due to the Lorentz force coming from the pseudo- 
magnetic field, the packet, which started in the K valley, 
is reflected by the upper armchair edge and scatters to 
the K' valley, where the pseudo-magnetic field points in 
the opposite direction. The semi-circular trajectory now 
travels in the opposite direction until the packet reaches 
the edge again, being scattered back to its former Dirac 
cone at the K valley. This procedure occurs several times 
until the packet is so strongly dispersed that it, eventu- 
ally, does not reach any of the edges, and performs only 
circular trajectories in the middle of the graphene flake. 
The series of reflections by the armchair border obtained 
in the strained case suggests the existence of a quasi- 
bound state at this edge, which is clearly seen by the 
time-dependence of the average coordinates (x) and (y) , 
shown in Fig. 5(b). As time elapses, both (x) and (y) 
simply oscillate around x = and close to the upper 
border of the sample, respectively. Notice that differ- 
ently from Fig. 4(a), the lateral and bottom borders of 
the panel in Fig. 5(a) do not match the edges of the sam- 
ple, in order to help the visualization of the trajectory, 
which in this case is localized in a small region of the sam- 
ple. One can also observe that the results for different 
wavepacket widths d are qualitatively the same, differing 
only by the amplitudes of the (x) and (y) oscillations in 
time. 
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FIG. 6: (Color online) Contour plots of the time evolution 
of the wavepacket in reciprocal space corresponding to the 
propagation shown in Fig. 5, i.e. for a strained graphene 
sample. The Fourier transform of the wavepacket is taken in 
the vicinity of the (a) K and (b) K' points of the reciprocal 
space illustrated in Fig. 3(a), along the k x — — 0.06A _1 (k x = 
0.045A -1 ) vertical axis for K (K'). 

The sequence of scatterings between K and K' val- 
leys suggested by the trajectories observed in Fig. 5 
is confirmed by a direct observation of the wavefunc- 



tion in reciprocal space. This is illustrated by Fig. 6, 
which shows the Fourier-transformed wavefunction (con- 
tour plots) as a function of the vertical component of the 
wave vector k y as function of time t, similar to Figs. 3 
(c) and (d). Since the initial wavepacket in this case is 
at k = (— 0.06 A _1 , 47r/3V / 3«), in the vicinity of the K 
point, the k y in Fig. 6(a) is taken for a fixed k x = —0.06 
A -1 , i.e. in the center of the initial wavepacket. One 
clearly sees that the peak of the wavepacket in reciprocal 
space oscillates between the K (a) and K' (b) regions in 
Fig. 6 as time elapses, as a consequence of the successive 
inter- valley scatterings by the upper armchair edge of the 
system, as observed in Fig. 5. Indeed, the first peak in 
the K' cone (Fig. 6(b)), for instance, starts to appear at 
t ~ 200 fs, which is the same time when (x) and (y) start 
to decrease in Fig. 5(b), after the wavepacket is scattered 
by the edge for the first time. It is important to point 
out that we had to take the k x = 0.045 A _1 vertical axis 
in order to observe the K' propagation in Fig. 6(b), in- 
stead of the k x = 0.06 A _1 that would be expected from 
the value of the wave vector in our initial wavepacket. 
In fact, one cannot expect that the K and K' points in 
the strained case remain vertically aligned in reciprocal 
space, as illustrated in Fig. 3(a), due to the distortion of 
the Dirac cones caused by the strain.^ Also, we observe 
that the scattered wavepackets are no longer gaussian, 
as they start to exhibit interference patterns due to the 
scattering by the edges. Even so, the conclusions drawn 
from the results in Fig. 6 are not affected by this fact, 
while Fig. 6 gives us a clear demonstration that the suc- 
cessive KtoK' scatterings are indeed strongly related to 
the non-propagating edge states found for the strained 
case in Fig. 5. 



IV. PROBABILITY DENSITY CURRENT 
CALCULATIONS 

The trajectories illustrated in Figs. 4 and 5 have a di- 
rect effect on the probability density currents, which are 
numerically calculated based on the method developed 
in Ref. [33j Since we can define the probability current 
j in terms of the continuity equation, then the discrete 
current centered on site n can be written as 

9 

Jn -Jn+l = CL-^Pn^n, (2) 

where p n ^ n = (n\p\n) are the matrix elements of the den- 
sity matrix operator p = and the time derivative 
is determined by the equation of motion for p 

d i 

di Pnn = li^ (*n*™-ffmn - H nm * m V* ) (3) 

m 

where \£ n = (n|\I/). We will limit ourselves to the case of 
nearest-neighbor interaction, i.e. i7 n ,™ = when \m — 
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| > 1, from which we obtain 
d i 



h 



[¥„¥;_!#„-!,„ - H„,n-l*n-l*;] , (4) 



which is easily rewritten in the form 



d 



3 [* n ** n+1 H n+1 , n ] 

9[* n _i*;fr n , n _i]. 



(5) 



By comparing Eqs. Q and (JEJ, one easily identifies the 
local current in n as 



2a 



3 



(6) 



Notice that Eq. ([6| was developed without taking into 
account any specific lattice and the presence of magnetic 
fields. However, a generalization to arbitrary discrete 
lattice is straightforward, and the presence of a magnetic 
field is included simply by the Peierls substitution of the 
hopping parameters.^ As graphene is a hexagonal lat- 
tice, the current components in x and y directions have 
different forms and are site dependent. Defining the sites 
location through their line (n) and column (m) positions 
in the lattice (see Ref. 29), one obtains 



j x (n,m) = ±- {23 [^n,m^n,m±l r n,m±l] 
-3 [*n l m*n-l 1 m^n-l,m] " 9 [*n,m* n+i >m Tn+l,m] } (?) 



and 



-3[*n,m*;_l, m r n _i |m ]} 



(8) 



where the =p sign in will be positive (negative) if the 
(n, ra)-site belongs to the sublattice A (B), and r n ^ m 
is the hopping parameter which, in the presence of a 
magnetic field, includes an additional phase according to 



the Peierls substitution r n , 



exp 



r a • di 

J rn 



where A is the vector potential describing the magnetic 
field. 

The probability density currents calculated by Eq. (|6| 
are integrated in space and plotted as a function of time 
in Fig. 7 for the situations proposed in Figs. 4 and 5, i.e. 

(a) in the presence of an external 5 T magnetic field, and 

(b) in a strained graphene sample, which produces an 
almost uniform « 5 T pseudo-magnetic field. As in the 
other results discussed previously, the results obtained 
for the two different values of wavepacket width consid- 
ered in this case, d = 100 A (black solid - circles) and 
d = 300 A (red dashed - triangles), exhibit similar qual- 
itative features, differing only in a quantitative way. In 
Fig. 7(a), for an external field, one observes a total cur- 
rent flow in the ^-direction oscillating around a positive 
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FIG. 7: (Color online) Integrated probability density currents 
as a function of time for the situations proposed in Figs. 4 
and 5, namely, (a) for an unstrained graphene sample in the 
presence of an external 5 T magnetic field, and (b) in a cir- 
cularly bent graphene sample, which produces an almost uni- 
form w 5 T pseudo-magnetic field. The curves (symbols) 
represent the component of the current in the x (^-direction, 
j x (j y )- Two different values of wavepacket width are 
considered: d = 100 A (black solid - circles) and 300 A (red 
dashed - triangles) The regions delimited in (a) are the same 
as in Fig. 4. 



value in the region I, whereas j y oscillates around zero in 
this region. This is a manifestation of the propagation 
of the wavepacket through the upper edge of the sam- 
ple, by means of skipping orbits, as illustrated in Fig. 
4(a). What follows can also be understood by analyz- 
ing Fig. 4(a): In regions II and IV (III and V), where 
the wavepacket propagates along the vertical (horizontal) 
edges, the component of the current in the 7/(x)-direction 
oscillates around a non-zero value, indicating an electron 
propagation through the sample by the skipping orbits 
mechanism. This is not the case when we consider a 
strain-induced pseudo-magnetic field: Fig. 7(b) shows 
that both j x and j y always oscillating around zero, con- 
firming that there is no net current in the system and 
that the skipping orbits near the armchair edge in this 
case are non-propagating states. 

The results found in our work have observable conse- 
quences in experiments. For example, the edge propaga- 
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tion of electrons through skipping orbits in an ordinary 
system under an external magnetic field plays an impor- 
tant role in electron transport in the direction parallel 
to the edge.^ Our results demonstrate that these skip- 
ping orbits are still present in a graphene ribbon under 
an external magnetic field, but they are not observed 
in the case of pseudo-magnetic fields in the direction 
parallel to an armchair border of graphene. This is a 
clear example that the pseudo-magnetic field has differ- 
ent consequences as compared to a real magnetic field. 
Therefore, in a strained armchair graphene ribbon, edge 
electrons should not propagate along the ribbon, so that 
the transport in these systems must be dominated only 
by electrons propagating far from the edge. Moreover, 
the non-propagating state found at the armchair edge 
of a strained sample is a consequence of periodic inter- 
valley scattering processes, and this type of scattering has 
an important effect on Raman spectroscopy. 35 Therefore, 
the successive electron reflections at the armchair edge 
of a strained sample would manifests itself as an intense 
peak in Raman experiments taken close to the border of 
the graphene sample. 

V. CONCLUSIONS 

In summary, we investigated the reflection of a 
wavepacket on zigzag and armchair edges of a graphene 
ribbon. Our results demonstrate the scattering of the 
wavepacket from K to K' Dirac cones in the case of arm- 
chair edges, whereas scattering from positive to negative 



average momentum inside the same cone is observed in 
the zigzag case, which is in agreement with predictions 
from mean field (Dirac) theory of graphene and with re- 
cent experimental results.^ In the presence of an external 
magnetic field, skipping orbits are observed. However, 
for a strain induced pseudo-magnetic field, our numeri- 
cal results demonstrate that the incoming and scattered 
wavepackets perform orbits in opposite directions in the 
armchair case. This effect is easy to be understood if one 
considers the combination between two events, both al- 
ready predicted by the Dirac theory of graphene: (i) the 
K to K' scattering by armchair edges and (ii) the opposite 
sign of the pseudo- magnetic field in the different cones. 
This result points directly to the possibility of observing 
non-propagating edge states in an armchair terminated 
strained graphene sample under pseudo-magnetic fields, 
which is completely different from the external magnetic 
fields case, where the skipping orbit states are always 
propagating. The effects predicted by our theoretical 
work are expected to have important consequences in fu- 
ture experiments on strained graphene samples. 
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